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^ ! ABSTRACT 

' We exploit the gravitational potential of massive cluster lenses to probe the emission 

\ line properties of six z — 1 galaxies which appear as highly magnified luminous arcs. 

. Using the GMOS integral field spectrograph together with detailed cluster lens models 

' we reconstruct the intrinsic morphologies and two-dimensional velocity fields in these 

galaxies on scales corresponds to ~0.5kpc (unlensed) at z = 1. Four of the galaxies 
^ ■ have stable disk-like kinematics, whilst the other two resemble interacting or starburst 

\->i \ galaxies. These galaxies lie close to the mean rest-frame /-band Tully-Fisher relation 

, for nearby spirals suggesting a clear preference for hierarchical growth of structure. In 

■ the rest-frame i?-band, the observations suggest 0.5±0.3mag of brightening, consistent 

^ \ with increased star-formation activity at z = 1. However, the galaxies with stable disk 

• 1^ . kinematics have more slowly rising rotation curves than expected from galaxies with 

/\ ' similar surface brightness in the local Universe. We suggest that this may arise because 

H I the distant galaxies have lower bulge masses than their local counter-parts. Whilst this 

. - - 1 study is based on only six galaxies, the gain in flux and in spatial resolution achieved 

via gravitational magnification provides a much more detailed view of the high redshift 
Universe than possible with conventional surveys. 

Key words: galaxies, Tully-Fisher relation, gravitational lensing, galaxy clusters. 
Integral Field Spectroscopy, Gravitational Arcs: Individual 



1 INTRODUCTION 



Massive galaxy clusters magnify the light from galaxies that 
serendipitously lie behind them. This natural magnification 
provides the opportunity to study intrinsically faint high 
redshift galaxies with a spatial resolution and to surface 
brightness limits that cannot be attained via conventional 
observations. These highly magnified sources can provide 
unique insights into the properties of typical galaxies at early 
times in the U niverse (Sniail ct al. 1996; Franx ot al. 1997; 
Teplitz et alTi oOO: Ellis ct al]|200ll:Tca mpusa no et alllioOlj 
Smith et al.ll2002. : .Swinbank et aljl2003i: iKneib et al.ll2004l) . 
The majority of known galaxy cluster lenses are at z ^ 0.3, 
which makes them ideal for detailed studies of galaxies at 



z ;> 1 jEdee et alJ 120031: ISand et alJl2005ll . i.e. before the 
Universe had reached half of its current age. 

At 2 ^ 1 important questions regarding the relation- 
ship between the total mass of a galaxy and the bary- 
onic mass locked up in stars remain unanswered. In lo- 
cal rotationally-supported spiral galaxies, this relation is 
best described by the Tully-Fisher relation iTuUv fc Fished 
Il977l) - an empirical correlation between the terminal ro- 
tational velocity (or line width) and the absolute magni- 
tude of a spiral galaxy. This relationship may reflect how 
rotationally-supported galaxies formed, perhaps suggesting 
the presence of self regulating processes for star formation in 
galactic disks. Measuring the evolution of the Tully-Fisher 
(TF) relation as a function of look back time therefore pro- 
vides important insights into the growth of galaxy mass 
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JVoet et al]ll997l: iBohm et aljlioo^ : iMilvang- Jensen et alJ 
l200a: iBamford et alJl2005h . One specific and rather simple 
test is to determine whether the steUar mass in galaxies 
has increased in lock-step with the mass of the dark matter 
halo or whether the stellar mass grows within a pre-existing 
dark matter halo. The former is pr edicted by "hierarchi- 
cal" galaxy formation models fe.e. [ White fc Frenld Il99ll : 
iKau ffmann e t aljUool ICole et alJl200d) since the growth 
of the stellar disk is regulated by the rate at which gas is 
accreted by the halo (the accretion rate is similar for gas and 
dark matter). The latter is characteris tic of simple ''classi- 
cal" galaxy formation models such as lEggen et al.l il962l) . 
In simplified outline, the hierarchical model predicts that 
the correlation between the terminal rotation velocity and 
stellar mass should evolve little, while the classical model 
predicts that the stellar mass corresponding to a fixed rota- 
tion velocity should decrease with increasing redshift. 

This observational test can be made using the galaxies 
photometry as a proxy for stellar mass. It is now generally 
accepted that the shape, zero-point and scatter of the TF 
relation depend on bandpass and that a smaller dispersion 
is obtained for near-infrared band-passes (i.e. rest-frame / 
rather than rest-frame _B-band). The dust-correction is much 
smaller in near-infrared photometry and moreover, it is more 
sensitive to the underlying, evolved stellar population which 
best traces the stellar mass and thus correlates more tightly 
with the maximum rotational speed (or total galaxy mass; 
IVerheiienl200ll : IConselice et alj2005l) (we note that, for gas- 
rich systems the bluer pass-bands are more sensitive to cur- 
rent star-formation and therefore may also correlate closely 
to total baryonic mass). 

However at z^l the small angular size of galaxies means 
that obtaining spatially resolved rotation curves is extremely 
challenging (e.g. a typical local spiral has a scale length of 
4kpc, which corresponds to only 0.5" at z—1). One solution 
to this problem is to exploit the magnifying power of fore- 
ground galaxy cluster lenses - for a typical lens magnifica- 
tion factor of ten, 0.6" corresponds to an unlensed physical 
scale of just ;£0.5kpc at z = 1. Thus galaxies can be tar- 
geted that would otherwise be too small. Spatially resolved 
kinematics can also be achieved on spatial scales far greater 
than otherwise possible due to the increased fiux sensitiv- 
ity over conventional observations. Moreover, the benefits 
of gravitational magnification are complemented by Inte- 
gral Field Spectroscopy (which produces a contiguous x, y, 
velocity map at each point in the galaxy) of the target galax- 
ies. Clean decoupling of the spatial and spectral information 
is therefore feasible, thus eliminating problems arising from 
mixing of the two in traditional long-slit observations. It is 
therefore much easier to identify which galaxies have regu- 
lar (bi-symmetric) velocity fields for comparison with local 
spirals. 

In this paper we present a study of six z=l 
gravitationally-lensed galaxies observed through the cores 
of four galaxy cluster lenses. Five of the targets lie be- 
hind Abell2390, CI 2236-04 and RGB 1745-^398, the sixth 
(Arc#289 at z=1.034 b ehind A 2218) was previously dis- 
cussed in ISwinbank et a l. ( 2003) . All six targets were ob- 
served with the Gemini Multi-Object Spectrograph In- 
tegral Field Unit (CMOS IFU) on Gemini-North. We 
concentrate on the galaxy dynamics as traced by the 
[Oll]AA3726.1,3728.8A emission line doublet. The IFU data 
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Table 1. Redshifts for the cluster lens (z^i); arcs {zarc) and spec- 
troscopic imaging exposure times for the arc s our sample. Ref er- 
ences for arc disco ve ry or previous study: l: |PeIIo et al.lll992l: 2: 
lEbbels et alll99^ 3 : ISwinbank et al.l2003l: 4:lNiIsson et al.ll999t 
5: iPello et al.l 119991: 6: iFrve fc Broadhurs j Il998t 7: iPello et alJ 
Il99lt 8: lKneib et al.lll993 



provide a map of the galaxy's velocity field in sky co- 
ordinates. To interpret this field we correct for the magnifi- 
cation and distortion caused by the lensing potential using 
models of the cluster lenses. These models are constrained 
by the positions and redshifts of other spectroscopically con- 
firmed gravitational arcs in each cluster. In a few cases, the 
IFU targets are themselves multiply imaged and the fold- 
ing of the velocity field places additional constraints on the 
cluster potential. These lensing corrections allow us to re- 
construct the intrinsic (unlensed) properties of the galaxies 
at z=l, including their geometry. 

The velocity field of systems displaying regular (bi- 
symmetric) rotational velocity fields i.e. resembling rotat- 
ing disks are reduced to traditional one-dimensional rota- 
tion curves. The terminal rotation velocity, indicative of the 
galaxy mass, is then extracted. All of the targets have been 
imaged in the optical and near-infrared bands. The grav- 
itational lens models allow us to correct the photometric 
observations for the amplification of the cluster and to es- 
tablish the intrinsic (unlensed) rest-frame B- and /-band 
luminosities of the galaxies. We can then compare the TF 
relations for our target galaxies with local observations, and 
hence to test the models for the build-up of galaxy mass as 
a function of look-back time. 

The structure of this paper is as follows. In JHwe de- 
scribe the sample selection, observations and the data re- 
duction. In we analyse the data, including photometry, 
gravitational lens modeling and galaxy reconstruction. In S^l 
we present the one-dimensional rotation curves and the TF 
relation at z=l. Finally in In f|3we discuss our results on 
the TF relation at z=l and the implications for galaxy evo- 
lution. We also outline the wider applicability of the grav- 
itational telescope method. Through-out this paper we use 
a cosmology with Ho = 72kms~^, Q,o = 0.3 and Ao = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample Selection 

The positions and redshifts of the gravitational arcs ob- 
served during this program are listed in Table 1. In order 
to avoid possible biases, the targets were selected to be rep- 
resentative of lensed galaxies in the distant Universe — no 
attempt was made to select galaxies with relaxed late-type 
morphology. However, we did require that arcs were resolved 
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Figure 1. Finding charts for the IFU observations. Top: True colour HST VI image of the lensing cluster A 2390 at 2=0.228 with the 
GMOS IFU field of view overlaid on the three 2~1 arcs. The yellow curve is the 2=0.912 critical curve calculated from the best-fit lens 
model of this cluster (' H3.2.2I see also Fig.|3J. We also label the z=4.04 and z=4.05 multiple-image systems that were used to constrain 
the lens model: H3-a/b/c and H5-a/b respectively. Lower left: True colour ground-based BVR image of the lensing cluster CI 2236-04 
with the GMOS IF U field of view overla id on the z=1.116 arc. The brightest two cluster galaxies are labelled A and B, and the 2=1.334 
source discussed bv iKneib et alJ 1 1994ft is labelled C. We also show the 2=1.116 critical curve in yellow. Lower right: True colour 
ground-based BRK image of the lensing cluster RGB 1745-1-398 with the GMOS IFU field of view overlaid on the 2=1.056 arc and the 
2=1.056 critical curve shown in yellow. The multiply-imaged ERO at (-5",-9") has a redshift of 2=1.11±0.05 ( i|2.4i and is labelled E\, 
E2 and E^. 
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in both spatial dimensions so that a two dimensional velocity 
field could be extracted from the IFU data. This restricted 
our selection to galaxies with moderate magnification. 



2.2 Imaging 

Near-infrared imaging of Abell 2390 was obtained th e Wide- 
field Infrared Camera (WIRC: IwTlson et alj liool on the 
Hale 200" Telescope^ . The J-band data were taken on 2004 
June 30, totaling 12.8 ks. The isT-band data were taken over 
four observing runs between October 2003 and June 2004 
yielding a total integration time of 11.5 ks. All of these 
data were taken when transparency was good, with the J- 
band conditions likely photometric. The seeing was typically 
% 0.85". The data were reduced using standard iraf^ proce- 
dures. Photometric calibration was achieved using 2MASS'^ 
data on stars in the same fields. We also retrieved observa- 
tions of A 2390 from the Hubble Space Telescope (HST) pub- 
lic archive''. The Wide Field Planetry Camera 2 (WFPC2) 
observations through the 7gi4- and Vsss-band filters totalled 
10.5 and 8.4 ks respectively - the data were reduced using 
the standard STSDAS package in iraf. The final reduced data 
are shown as a colour image in (Fig.0. A 2390 has also been 
observed with the Advanced Camera for Surveys (ACS) as 
part of the Guaranteed Time Observations (GO: 9292, PI: 
Ford). We concentrated solely on the z'-band ACS data, re- 
ducing the single orbit (2.2 ks), CR-split, undithered data 
using MULTIDRIZZLE with the default parameter set. 

Ground-based optical imaging of CI 2236-0 4 and 
RG B 1745-F398 (Fig. [Tl is taken from iKneib et al.) lll994l) 
and lNilsson et al] |S93) respectively. These consist of B, R 
and _B, _R, 7-band imaging of these clusters taken with 
the New Technology Telescope (NTT) and Nordic Optical 
Telescope (NOT) respectively. We supplement these opti- 
cal data with near-infrared J- and A'-band imaging of both 
clusters on the UK Infra-Red Telescope (UKIRT) between 
2004 August 15 and 2004 August 20^. The observations were 
made in photometric conditions and ~ 0.7" seeing using 
the UKIRT Imaging Spectr ometer (UIST) imaging camera 
l|R,amsav Howat et al.Ml998^ which employs a 1024 x 1024 
InSb detector at 0.12"pixel~^ to give a 2' field of view. 
The observations were taken in a standard nine-point dither 
pattern and reduced using the relevant ORAC-dr pipeline 
iCavanagh et al.ll2003ll . The total integration times for each 



^ The Hale Telescope at Palomar Observatory is owned and op- 
erated by the California Institute of Technology. 
^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under a cooperative agreement with 
the National Science Foundation. 

This paper makes use of data products from the Two Micron 
All Sky Survey (2MASS), which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation. 
* Obtained from the Multimission Archive at the Space Telescope 
Science Institute (MAST). STScI is operated by the Association 
of Universities for Research in Astronomy, Inc., under NASA con- 
tract NAS5-26555. Support for MAST for non-HST data is pro- 
vided by the NASA Office of Space Science via grant NAGS- 7584 
and by other grants and contracts. 



band was 2.4 ks. To calibrate our data, we observed UKIRT 
faint photometric standards llHawarden et al.]l200ll) . These 
standards were observed at similar air masses and using the 
same instrumental configuration as the target galaxies. 

The most striking feature in the new near-infrared imag- 
ing is the discovery of a red triply-imaged galaxy 4" to the 
South-East of the z=1.056 arc in RGB 1745-^398 (Fig. Ql. 
This arc is barely detected in the optical imaging, but is 
very bright in the near-infrared (J = 18.27 ± 0.06 and 
K = 16.32 ± 0.04) and &nR-K colour of 6.6 ± 0.2. The ob- 
vious triple image configuration, along with a spectroscopic 
redshift of this arc would provide a new constraint on the 
lens model of this cluster and therefore a more precise recon- 
struction of the intrinsic properties of the blue arc targeted 
by the IFU observations ( H2.3|l . We describe ground-based 
optical and near-infrared spectroscopy of the red arc in H2.4I 

2.3 Integral Field Spectroscopy 

Traditionally, observations are made with a two-dimensional 
detector, sufficient for imaging programs. However, for spec- 
troscopic measurements, one spatial dimension of the sky is 
usually lost in order to allow dispersion of the image across 
one axis of the detector. Spectrograph's placing a long slit 
over the object lose spatial information orthogonal to the 
slit width. In order to faithfully investigate the properties 
of galaxies which are spatially extended, we require spec- 
troscopic information over the whole of a two-dimensional 
field. An instrument simultaneously producing both two- 
dimensional imaging data and a spectrum for each point 
in the image provides integral field spectroscopy. 

Integral field spectroscopy thus allows spectra to be 
simultaneously obtained for a number of contiguous areas 
across a two-dimensional field. The spatially resolved spec- 
troscopy allows the variation with positions of object spec- 
tra to be measured, and this information can be used, for 
example, to determine the object's internal kinematics, or 
the star formation rate within an object, providing insight 
into its process of evolution. This technique is therefore cru- 
cial in order to investigate the internal properties of distant 
galaxies. 

Our sample were observed in queue mode with the Gem- 
ini Muhi-Object Spectrograph (CMOS) on Gemini-North'' 
between 2003 June 8 and 2003 August 14 in photometric 
conditions and gO.7" seeing. Details of the observations and 
source redshifts are listed in Table 1. The IFU uses a lensed 
fiber system to reformat the a contiguous 7.9" x 5.3" field 
(comprising an array of 40x2 6 hexagonal fibres each of 0.2 " 
diameter) into two long slits dAllineton-Smith et aljl2002ll . 
All observations were made using an /-band filter in conjunc- 
tion with the R400 grating which results in two tiers of spec- 
tra covering the maximum field of view. The spectral resolu- 
tion of this configuration is A/AA = 2000. Each observation 
was split into 2.4 ks sub-exposures and dithered by one IFU 
lenslet to account for bad pixels. The data was reduced in 
IRAF using the GMOS-IFU data reduction pipeline which 
extracts, flat-fields and wavelength calibrates the data. To 

® The United Kingdom Infrared Telescope is operated by the 
Joint Astronomy Centre on behalf of the U.K. Particle Physics 
and Astronomy Research Council. 
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Figure 2. Top: FST/WFPC2 and CMOS IFU observations of A2218 arc#289. Left: HST/WFPC2 B450V6O6/814 image of A2218 
arc#:289. Middle: The [On] A3727 emission map of the arc measured from our IFU observations. The distribution of [On] emission agrees 
well with the rest-UV flux seen in the left panel. Right: The velocity field of the galaxy derived from the [On] emission. The scale is 
marked in arcseconds and North is up and East is left. Bottom: The reconstructed galaxy after correcting for lens magnification — see 
H3.2.1l for details. Left: the reconstructed image of the galaxy. Middle: The reconstructed [On] emission line map. Right: The reconstructed 
velocity map. The red and blue regions represent redshift and blueshift respectively, and the solid lines show the asymptotic major axis 
cross section from which the one-dimensional rotation curve was extracted (' H4.1i . The scale shows the size of the galaxy after correction 
for lensing. 



subtract the sky emission lines the GMOS IFU employs a 
second (5" x 3") IFU separated by one arcminute on the sky. 
We used idl to identify and extract sky-fibres adjacent to ob- 
ject fibres on the spectrograph and used these to achieve the 
sky-subtraction. We also improved the flattening of the data 
by using continuum regions either side of the emission lines. 
This is achieved by identifying wavelength regions free from 
sky emission lines near the redshifted [On] line and averaging 
the flux in these regions over 200A. For each fiber, this aver- 
age was divided into the fiber spectrum to improve the fiber- 
to-fiber response map. The output pixel scale is I.SA pixel 
and the instrumental profile has a FWHM of 3.4A (mea- 
sured from the widths of the skylines). This corresponds 
to a = 58kms~^ in the galaxy rest-frame, and the set-up 
adequately resolves the [Oii]AA3726.1,3728.8A emission-line 
doublet. In all following sections we have deconvolved the 
instrumental resolution from the line width measurements 
where these could be reliably determined. 



Programme ID: GN-2003A-Q-3. The GMOS observations are 
based on observations obtained at the Gemini Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on 



For one of the targets, (A2390arcA), we required two 
pointing's to cover the large spatial extent (~ 12" x 4") of 
the galaxy. During the second pointing the much smaller arc 
(~ 2" X 2" at z =1.033; A 2390arcB) was also covered. In or- 
der to ahgn and mosaic these two datacubes we constructed 
(wavelength-collapsed) white-light [On] emission line maps 
of the two arcs. This was achieved by identifying the central 
wavelength of the [On] emission in each galaxy and then 
collapsing the data-cubes between -300 and -1-300 km s~^. 
These were then compared and aligned using the HST I- 
band image as a reference to produce the final mosaiced 
datacube. 

We also note that the redshift of the arc in 
RGB 1745-1-398 (z=1.056) places the [On] emission at 
7660A; close to the Fraunhofer A-band. To correct for this 
absorption, we extracted the spectrum of the foreground 
cluster galaxy (which is a strong continuum source at this 



behalf of the Gemini partnership: the National Science Founda- 
tion (United States), the Particle Physics and Astronomy Re- 
search Council (United Kingdom), the National Research Council 
(Canada), CONIC YT (Chile), the Australian Research Council 
(AustraUa), CNPq (Brazil) and CONICET (Argentina). 
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Figure 3. Top left: HST and IFU observations of A2390arcA and arcB. We show tlie V606^814- band image of tlie two arcs (left), 
tile [On] emission map (middle) and the velocity field of both galaxies (right). The distribution of [On] emission in both galaxies agrees 
well with the rest-UV flux seen in the HST imaging. In all three panels the 2=0.912 critical curve is shown in yellow/black. Top right: 
The reconstructed galaxies after correcting for lens magnification - see i|3.2.2l for details. The upper and lower three panels show the 
reconstructed image, [On] emission line map and the velocity field of arcA and arcB respectively. In the velocity maps the red and blue 
regions represent redshift and blueshift respectively, and the solid lines show the asymptotic major axis cross section from which the 
one-dimensional rotation curve was extracted. Middle: HST and IFU observations of A2390arcD. Far left: V606^8l4-band color picture 
showing the complex morphology of this red galaxy. Left center: The [On] emission line map. The distribution of [On] emission line flux is 
well matched to the broad-band imaging, and confirms that the bright knot (located 3" to the south) is associated with this galaxy. Right 
center: The velocity field derived from the [On] emission. Far right: The [On] emission line width (FWHM) map which shows a large 
variation in the distribution of line widths (both in the spatial and spectral domain). Bottom: The reconstructed galaxy after correction 
for lens magnification. The galaxy has a complex morphology, with at least two components. The peak-to-peak velocity gradient of the 
galaxy is orthogonal to the major axis seen in the broad-band imaging. Furthermore, the distribution of the [On] line widths show two 
(off-center) peaks approximately along the same direction as the maximum velocity gradient. The dense knot or companion located 7 
kpc to the south appears redshifted by 300±80kms~^ and has a line width of ~300kms~^ FWHM. 
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Figure 4. Top: Ground-based imaging and IFU observations of the 2=1.116 arc in CI 2236-04. Left: BV R-ha,nd colour picture. Middle: 
[On] emission line map . Right: The velocity field derived from the [On] emission. In all three panels the 2=1.116 critical curve of 
the best-fit lens model is shown in yellow/black. Bottom: The reconstructed galaxy after correction for lens magnification. Left: the 
reconstructed BV R-ha,nd colour picture. Middle: The reconstructed [Oil] emission line map. Right: The reconstructed velocity map. The 
red and blue regions represent redshift and blueshift respectively, and the solid lines show the asymptotic major axis cross section from 
which the one-dimensional rotation curve was extracted. 
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Figure 5. Top: Ground-based imaging and IFU observations of the 2=1.056 arc in RGB1745+398. Left: Ground-based colour BRK- 
band colour picture. Middle: [Oil] emission line map. Right: The velocity field derived from the [On] emission. In all three panels 
the 2=1.056 critical curve from the best-fit lens model is shown in yellow/black. Bottom: The reconstructed galaxy after removal of 
lens magnification. Left: Reconstructed BRK-hand colour picture. Middle: The reconstructed [On] emission line map: note that the 
morphology seen in this image is largely determined by the IFU field of view. Right: The reconstructed velocity map. The red and blue 
regions represent redshift and blueshift respectively, and the solid lines show the asymptotic major axis cross section from which the 
one-dimensional rotation curve was extracted. 



wavelength; see tl3.2.4ll and use this to model and correct for 
the telluric absorption. 

After constructing the datacubes, we proceed to fit the 
[On] emission line doublet in each pixel using a min- 
imisation procedure, taking into account the greater noise 
at the position of the sky lines. The spectra were averaged 



over 3x3 spatial pixels (0.6" x 0.6"), increasing this region 
to 4x4 pixels (0.8" x 0.8") if the signal was too low to give a 
sufficiently high improvement over a fit without the line. 
In regions where this averaging process still failed to give 
an adequate i no fit was made. Using a continuum fit, we 
required a minimum oi 25 (S/N of 5) to detect the line 
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Figure 6. LRIS-R spectra of the multiply imaged ERO in 
RGB 1745+398. The spectrum is produced by rejecting pixels 
dominated by the strong night sky emission i.e. those regions 
above the horizontal line in the lower spectrum which shows the 
night sky spectrum (offset and scaled for clarity). The smoothed 
spectrum is produced by median filtering the raw spectrum with 
a boxcar width of lOOA and the error bars represent Itr errors in 
each boxcar window. We identify the discontinuity in the spec- 
tral shape at A8440A as the 4000A break which yields a redshift 
of 2=1.11 ±0.05. The blue line shows an Luminous Red Galaxy 
(LRG) galaxy template from the Sloan Digital Sky Survey red- 
shifted to 2=1.11. 

and when this criterion is met, we fit the [On] emission line 
doublet with a double Gaussian profile of fixed separation, 
allowing the normalization and central wavelength to vary. 
If the line is detected at greater than 7a significance, we also 
allow the line width to vary. The two lines of the doublet 
are assumed to have the same intensity and width. Whilst 
it is possible for the two emission lines of the [On] doublet 
to have different intensities (due to the density and tem- 
perature of the gas), at this spectral resolution and typical 
signal-to-noise we view it as more reliable to fix the ratio 
of the intensities to unity (we note that in the collapsed 
spectrum from each galaxy there is no noticable difference 
between in the fluxes between the two peaks). To calculate 
the error in the velocity we perturb the wavelength of the 
best fit profile and allow the signal to drop by a of 9. 
This corresponds to a formal Scr error. 

2.4 Long Slit Spectroscopy 

In the course of assembling the new and archival imaging 
data required for this article we discovered a new triply- 
imaged red arc in RGB 1745-f398 (Fig.0. We estimate the 



redshift of this arc using hyper-z iBolzonella et al.ll200oll 
and the photometry given in H2.2I obtaining z = 1.1 ± 0.2. 
We attempted to measure a spectroscopic redshift for this 
arc on 2005 February 22. We used the Cooled Grating Spec- 
trograph 4 (CGS4) and UIST, both on UKIRT to explore 
1.0— 1.2/im and 1.4— 2.4/xm respectively - the total integra- 
tion times for both of these observations was 1 hour. We de- 
tected continuum emission in these data, but no obvious line 
emission. We therefore also searched for spectral features at 
shorter wavelengths using the r ed arm of the Lo w Resolution 
Imaging Spectrograph (LRIS: loke et alJll995h in long slit 
mode on the Keck I 10-m telescope'^ on 2005 March 05. The 
831/8200 grating was oriented at 36° to sample the wave- 
length range ~0.75— 0.9/im. A total integration time of 3.3 ks 
was accumulated from three separate exposures with a 1" 
wide slit. Flux calibration was achieved via observations of 
the spectrophotometric standard G138-31 (Oke 1990). The 
data were reduced using standard iraf tasks - the final re- 
duced 1-dimensional spectrum is shown in Fig.|S| Again, no 
strong emission or absorption features were found, however 
there is a continuum break at 8440A which we interpret as 
the 4000A break redshifted to z = 1.11 ± 0.05. The uncer- 
tainty on the redshift is estimated by comparing the data 
with the mean spectrum of a luminous red galaxy from the 
Sloan Digital Sky Survey; (Fig. [SJ. 



3 DATA ANALYSIS AND MODELLING 
3.1 Photometry 

From our optical/near- infrared imaging, we constrain the 
spectral energy distribution (SED) of each galaxy. Since the 
arcs usually lie with a few arcs-seconds of nearby bright clus- 
ter galaxies, we calculate the magnitude of the arcs in var- 
ious pass-bands by masking the arc and interpolating the 
light from the nearby cluster members. We then use the 
IRAF IMSURFIT package with the sky estimated from a 2"^* 
order polynomial surface fit. This is subtracted from the 
regions around each galaxy to account for the halo light 
from the nearby cluster galaxies. We then use SEXTRAC- 
TOR (Bertin & Arnouts 1996^ to estimate the residual back- 
ground within the frame and extract the arc photometry us- 
ing SEXTRACTOR. We use the observed colours (at the known 
redshift) to infer the ratio of current to past star formation 
rate and use this to find the best fit SED. At 2: ~ 1 the rest 
frame /-band luminosity is approximately equivalent to the 
observed _ff-band and so is calculated from interpolating J 
and K photometry. We co-add the J, K images to obtain an 
aperture which is used to extract the J + K magnitudes and 
interpolate for the relevant SED type. We apply the same 
technique for the rest-frame B-band magnitude, (at z ~ 1, 
R and /-bands are the closest match to rest-frame /3-band) 
and estimate its uncertainty by computing the magnitude 
for a variety of SED types which are consistent with the 
observed optical colours. 

^ Obtained at the W.M. Keck Observatory, which is operated as 
a scientific partnership among the California Institute of Tech- 
nology, the University of California and the National Aeronautics 
and Space Administration. The Observatory was made possible 
by the generous financial support of the W.M. Keck Foundation. 
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Table 2: Arc Photometry 



A2218arc#289 21.66 ± 0.04 20.86 ± 0.10 20.53 ± 0.04 19.79 ± 0.05 18.63 ± 0.06 17.32 ±0.06 

RGB1745 + 398arc 21.90 ± 0.12 21.26 ± 0.10 20.90 ± 0.07 19.6 ±0.3 18.29 ± 0.09 16.72 ±0.07 

A2390arcA ... 21.64 ± 0.06 ... 19.91 ± 0.04 18.64 ±0.10 17.42 ± 0.08 

A2390arcB ... 23.34 ± 0.07 ... 21.65 ± 0.05 20.10 ±0.10 19.41 ± 0.06 

A2390arcD 24.55 ± 0.17 21.62 ± 0.04 ... 19.53 ±0.03 17.31 ± 0.03 15.64 ± 0.02 

C12236-04arc 20.50 ± 0.11 ... 20.05 ±0.11 ... 19.47 ±0.14 18.09 ±0.12 



Table 2. Optical and near- infrared photometry for the arcs in our sample. 



Table Source-Frame Properties of the Arcs 



Amplification Inclination v^^j Mj Disk Scale Lengths 

(fj.) (i) kms~l hy hj '^[OI 



A 2218 arc#2S9 
RGB1745 + 398arc 


4 g2 + 0.20 


64 ± 6 
60 ± 7 


206±17 
103±25 


-21.10±0.20 
-19.55±0.60 


-22.34±0.21 
-20.22±0.55 


2.4 


2.1 


2.2 


A2390arcA 




69 ± 4 


170±28 


-20.09±0.24 


-21.37±0.25 


2.0 


2.2 


2.6 


A2390arcB 


18 + 1 






-18.24±0.30 


-18.70±0.20 


0.8 


0.8 


1.1 


A2390arcD 








-21.77±0.20 


-23.27±0.15 








C12236'04arc 


,,, + 2.8 


45 ± 10 


98±30 


-19.79±0.22 


-20.51±0.23 


1.6 


1.3 


1.9 



Table 3. Notes: fi denotes luminosity weighted magnification. We estimate the l—rr uncertainty in fi by perturbing the parameters of 
the best fit lens model such that Ax^ = 1 and for each model we recompute the magnification of the arc. The error corresponds to the 
largest variation in magnification in the source plane photometry. The inclination of the galaxies are measured in the source-frame and 
the associated error-bars are calculated by determining the inclination in each pass-band and for each of the lens models described above. 
The rotation velocities (vmi) are corrected for inclination effects. The disk scale lengths are given for the reconstructed / and V-band 
images as well as the distribution of [On] fiux. 



3.2 Gravitational Lens Modelling and Galaxy 
Reconstruction 

In order to investigate the source-plane properties of tiie 
sources we must first correct for tiie distortion and magnifi- 
cation of the galaxy image by the cluster lens. Although this 
does not affect the sign or amplitude of the velocities, we 
must determine accurately the gravitational magnification 
of each galaxy so that we can recover its intrinsic geome- 
try and hence estimate the inclination angle. In this section 
we describe briefiy the construction of the cluster lens mod- 
els and their application to derive the intrinsic properties of 
the galaxies at z=l. T he lens model l ing te chniques are those 
describe d in de t ail by ISmith et alJ (|2TO^, orig inally devel- 
oped by iKneib' fm3), and further refined bv lKneib et alJ 
il99(il and Smith (2002). The primary constraints on the 
lens models are the positions and redshifts of spectroscopi- 
cally confirmed gravitational arcs in each cluster. The IFU 
data ( t|2.31 provide additional constraints because the com- 
ponents of the observed arcs that are identified as being 
multiple-images of the same region of the respective lensed 
galaxies must, within the uncertainties, have the same ve- 
locity and [On] line strength, in addition to the usual broad- 
band fiux and colour constraints. 



3.2.1 A 2218 

Arc #289 in A221 8 was discussed in detail by 
ISwinbank et alJ (|200^. Here we summarize the key 
details and use this arc to explain the methods applied to 
the other arcs. 

A 2218 is one of the best constrained strong lensing clus- 
ters, including three multiple-image systems that have been 
spectroscopically confirmed via their line emission. This 



cluster has been modeled on numerous occasions; Kneib et 
al. (1993; 1996; 2004); Ellis et al. (2001); Smith et al. (2005). 
We use the most recent of these to calculate the luminos- 
ity weighted magnification (/x). This was achieved by ray- 
tracing between the image- and source-planes to build up a 
map across the arc of how the observed flux relates to the in- 
trinsic (unlensed) fiux. The statistical error on /i was derived 
from the family of lens models which adequately reproduce 
the multiply-imaged arcs. For each model we then recom- 
pute the magnification of each arc. The error corresponds to 
the largest variation in magnification that we found via this 
method. The luminosity weighted amplification for this arc 
is ^=4.92lg jg, which translates to a boost in magnitude of 
Am = 1.7 ±0.1. 

The mapping between image- and source-plane co- 
ordinates described above was also used to reconstruct the 
intrinsic morphology of Arc#289 from the fiux map of the 
arc after lensing correction, (see Fig. 3). The reconstruc- 
tion reveals that Arc#289 is a blue disk-galaxy with much 
internal structure, resembling a late- type galaxy. This is es- 
pecially prominent in the B-band which samples the rest- 
frame UV and is therefore dominated by the star-forming 
Hll regions (Fig. 3). To measure the geometry of the disk 
we fitted ellipses to an isophote in the reconstructed fiux 
maps using the idl GAUSS2dfit routine (and assume an in- 
trinsically circular disk). The ellipticity is then e — 1 — b/a 
(where a and b are the major and minor axis of the ellipse) 
and the inclination, i, is cosi — b/a. The average axis ratios 
of the ellipses from the various passband are translated into 
inclinations and the error bars on i are computed by apply- 
ing the same procedure to each passband and each of the 
lens models described above. The inclination of the galaxy 
is found to be 64 ± 6° . 

We use the intrinsic velocity field (constructed in a man- 
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ner similar to that applied to the broad-band imaging above) 
to infer the rotational velocity of the galaxy's disk. The 
terminal rotation velocity is 186 ± 16kms~^ (i.e. half the 
asymptotic velocity shift across the galaxy). Using the in- 
clination angle derived above, we therefore calculate a cor- 
rected rotation velocity of Wrot = (206 ± 17) kms~^. 

3.2.2 A 2390 

A cluster lens mo del of A 2390 was originally developed by 
IPello et alJ il999|) who spectroscopically identified two z~4 
multiply imaged galaxies behind the cluster. The first (la- 
belled H3 in Fig. has a redshift of z = 4.04 and all three 
images lie within 10" of the cluster galaxy to the West. The 
second lies at z = 4.05 (labelled H5-a and H5-b in FigQ. 
However, the thir d image of this galaxy was not identified by 
iPello et alJ lll999ll because it lies outside of the field of view 
of the HST /WFPC2 observations available at that time. 
We have used the larger field of view of the new ACS z'- 
band data ( ^'2.21 to identify the third counter-image of H5 at 
Q=21:53:36.64, 5^-1-17:42:11.8 (J2000). These spectroscopi- 
cally confirmed multiple-image systems are used to constrain 
the lens model of this cluster, including the generation of 
a family of acceptable models as described in il3.2.1l The 
z = 0.912 critical curve of the best-fit model is shown in 
Fig. Q - this curve defines where multiply-imaged galaxies 
at 2=0.912 will appear in the observed (image) plane. In a 
simple axi-symmetric system the gravitational magnification 
is infinite along critical curves. This is not the case in the 
perturbed gravitational potential of a galaxy cluster, never- 
theless magnifications of ~10— 20x are routine (Table 2). 

A2390 arcA - This arc (the so-called "straight arc" - 
iPello et alJll999^ is -10" long, has a redshift oi z = 0.912 
and lies approximately 40" to the West of brightest clus- 
ter galaxy (BCG). Based solely on the broad-band imaging, 
the arc may be multiply-imaged at its southern end, in- 
deed the 2=0.912 critical curve of the best-fit lens model 
is immediately adjacent to the southern portion of this arc 
(Figs.0& 3). However the acceptable models include several 
that predict this arc to be strongly-sheared at its southern 
end rather than multiply-imaged. The velocity field derived 
from the CMOS observations supports the multiple-image 
interpretation: the southern end of the galaxy (So) has an 
observed velocity (with respect to the center of the galaxy) 
of ~ 150 km s~^ (Fig. 3) . The relative velocity of components 
5*1 and 5*2 also agree within la with the velocity of So. We 
conservatively include both multiple-image and strong-shear 
interpretations when estimating uncertainties below. 

We use the same methods as described in H3.2.1l to re- 
construct the intrinsic properties of this arc, obtaining the 
lensing corrected morphology, [On] map and velocity map 
shown in Fig. |31 This galaxy therefore appears to have disk- 
like kinematics, to be magnified by /i=12.6lQg, (Am = 
2.75 ± 0.07), and have an inclination angle of i=(69±4)°. 
The terminal velocity (uncorrected for inclination effects) is 
175 ± 20kms~^, which translates into a corrected rotation 
velocity of i;rot=(187±17) kms^"'. 

A2390arcB - The CMOS observations of A2390arcA 
also allowed us to study the [On] emission from a slightly 
higher redshift (z=1.033) arc, located 2" to the South 
of A2390arcA (Fig. [3 & OJ. This galaxy has a com- 



pact observed morphology (FWHMgl.5"). We calculate a 
luminosity-weighted magnification of 18.0 ± 1.0, (Am = 
3.1 ± 0.1) which yields a lensing corrected half-light ra- 
dius of just 2.5 kpc. There is no significant velocity gradient 
along the major axis of the reconstructed galaxy (Fig. 3), 
however we extract the maximum velocity shear from the 
data to derive an estimated a possible peak-to-peak ve- 
locity of 60±15kms~'^. Moreover, the [On] emission dou- 
blet is well resolved with a FWHM of ;£60kms~^ (decon- 
volved for instrumental resolution). These estimates place 
A2390arcB comfortably within the scatter of Hll regions 
in the velocity-width versus half light radius from lower 
redshift (z~0.19— 0.3 5) narrow-emission line galaxies from 
iGuzman et all Jl99el) . 

A2390arcD - This 2=0.912 arc lies 19" away from 
A2390arcA, their similar redshift suggesting that A 2390 
may be lensing a group of galaxies. The HST observations 
reveal an irregular morphology and the colours are consis- 
tent with either an evolved stellar population or a dusty- 
reddened star-burst. The strong and spatially resolved [On] 
emission revealed by our GMOS data appear to favor the lat- 
ter interpretation, and reveal that the galaxy is not relaxed 
(Fig. |3J. The mean velocity of the line emission suggests 
possible rotation about the major axis of the galaxy with 
peak-to-peak velocity of — 500 km s"'^. This is orthogonal 
to the direction of rotation suggested by the galaxy's mor- 
phology. Furthermore, the line profiles in the central region 
appear clearly broadened, having FWHM~300 kms~^. This 
suggests either (i) a merger between two galaxies, or (ii) 
shock heated, outflowing material from the central regions 
of the galaxy (powered by AGN activity). Our observations 
also confirm that the (unresolved) bright knot (located 3" 
to the south) is part of the same system, but show that it is 
offset by ~480±60kms~^ in velocity and has an [On] line 
width of 300±100kms"\ The magnification of this source 
is ^i=6.7l|5:2 = 2.0 ± 0.1). 



3.2.3 CI 2236-04 

The blue z=1.116 arc in CI 2236-04 jMelnick et alJll^^3^ 
lies between two bright elliptical galaxies (labelled A and 
B in Fig. 0. The ar c is almost straight with a length of 
8", and was shown bv lKneib et al.l il994f) to not be part of 
a multiple- image system including the blue sou rce labelled 
C in Fig. Q which is at 2=1.334. iKneib et al.l also identi- 
fied a velocity gradient along the 2=1.116 arc using longslit 
spectroscopic data; but the lack of full two-dimensional data 
lead them to suggest that the arc comprises a pair of closely 
interacting galaxies. The lack of both high resolution HST 
imaging and integral field spectroscopy of the arc made it 
difficult for them to test the galaxy-galaxy merger hypothe- 
sis. With the full two-dimensional coverage our GMOS data 
are therefore a powerful tool in this regard because they al- 
low us to decouple cleanly the spatial and spectral informa- 
tion that are mixed in long slit observations. Fig.^Jshows the 
[On] emission line intensity and velocity field of the arc. The 
velocity field across the arc is seen to be smooth and regular, 
with the eastern (receding) portion stretched by the grav- 
itational potential of the cluster. There is no obvious sign 
of a galaxy-galaxy merger in these data. The simplest inter- 
pretation is therefore that the arc comprises a single disk- 
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like galaxy, the receding portion of which is either multiply- 
imaged or strongly-sheared by the gravitational potential of 
the foreground cluster - i.e. a similar interpretation to that 
of arcA in A 2390. 

We construct a suite of nine models that are able to 
reproduce the observed arc morphology and velocity field. 
All of the models contain mass components for galaxies A 
(the BCG) and B, plus a mass component for the central 
cluster dark matter halo centered on the BCG. Seven of the 
models have bi-modal dark matter distribution and contain 
a dark matter halo centered on galaxy B. The models with a 
unimodal dark matter distribution (e.g. (TDAf,A=930 kms~^, 
o"DM, B=Okms~^) shear the galaxy at 2=1.116 sufficiently 
to reproduce a symmetric source-plane morphology with 
respect to the position where the observed velocity field 
changes sign. The more extreme of the bi-modal mod- 
els (e.g. (ToAf,A=750 kms~^, (jDM,s=500kms~^) multiply- 
image the receding part of the galaxy, but do not have quite 
enough mass to multiply-image the approaching part of 
the galaxy. The best fit lens model has (t_dj\/,a=875 kms~^, 
o"i3M,s=300 kms~^. The 2=1.116 critical curve of the best- 
fit model is shown in Figs. & |1] All nine models are used 
when estimating the uncertainties in the galaxy reconstruc- 
tion. 

Following the methods described in tl3.2.1l we use the 
suite of models to obtain a mean luminosity weighted mag- 
nification of fj.=11.2'^2'^ (Am=2.6±0.2) and a lensing cor- 
rected inchnation of i=(45±10)°. The observed velocity pro- 
file shows a terminal velocity of 68±15kms~^ across the 
galaxy which translates into an inclination corrected rota- 
tion velocity of Vrot=d8t.2(i) kms~^. 

3.2.4 RGB 1745+398 

RGB 1745-1-398 is qualitatively similar to CI 2236-04 as a 
gravitational lens in that a straight blue arc lies between the 
BCG and the second bright est cluster galaxy (F ig. [Jl. The 
arc was first discovered by iNilsson et al.l ^1999^ and spec- 
troscopically confirmed to be at 2=1.056. Our new CMOS 
IFU observations confirm that the bright source 1" to the 
North- West of the arc is a foreground cluster galaxy and is 
not part of the arc. 

We construct a suite of lens models to explore the 
balance of mass between the dark matter halos centered 
on galaxies A and B in a manner analogous to that de- 
scribed in il3.2.3l The model parameters are constrained 
by the triply-imaged red arc discussed in t|2.4l and the 
requirement to reproduce the observed morphology and 
velocity field of the blue arc. The best-fit model has 
o"DM,A=560 kms~^, (TDA/,s=400kms~^, although formally 
a very wide range of parameter space is allowed; the ex- 
trema models have (t_djv/,a=660 km s~^, (TDAf,B=90kms~^ 
and (TDA/,A=0kms~^, crDA/,s=630kms'^ respectively. We 
note that cr_DM,A=0 kms~^ is somewhat unphysical, how- 
ever rather than invoke a prior we conservatively use the 
full range of models allowed by the data alone. The critical 
curve of the best-fit model is shown in Fig. Q 

The velocity map of this galaxy exhibits a strong ve- 
locity gradient. The IFU only covered the central regions 
of the galaxy (as can be seen in Fig. with some low sur- 
face brightness emission lying outside the field of view. How- 
ever, the velocity field (and rotation curve) is smooth and 



characteristic of a rotating system and gives a good indica- 
tion of the terminal rotation velocity, which we determine 
to be 90±20kms~^. Using the family of models allowed by 
the data we derive a luminosity weighted magnification of 
/i=14.8lg g (Am=2.9lQ g), a lensing corrected inclination 
angle of z=(69±4)° and an inclination corrected rotation 
velocity of Vrot=W3'^23 kms^^. 



4 RESULTS 

4.1 One-Dimensional Rotation Curves 

Having reconstructed the source frame morphologies and 
velocity fields of the lensed galaxies, we can extract the 
one-dimensional rotation curves in order to determine their 
asymptotic rotation speeds as well as more general issues 
such as how the rotation curve shapes compare with those 
in the local Universe. Our sample consists of six galaxies, the 
fi ve presented in this pap er, plus A 2218 arc#289 presented 
in ISwinbank et al] i2003l) . In four of the lensed galaxies in 
our sample, A2390arcA, A 2218 arc#289, C12236-04arc and 
RGB 1745-|-398arc we identify and extract the major axis 
velocity profiles to measure the asymptotic rotation speed. 
This direction coincides with the maximum velocity gradi- 
ent as we would expect in a rotating system. The other two 
galaxies, A 2390arcB and A 2390arcD have velocity fields 
which are not aligned with the major axis. Although the ve- 
locity field is not consistent with a simple rotating disk, we 
are still able to place limits on their masses using the veloc- 
ity offsets between merging components (A2390arcD) and 
the limits on velocity gradients/line widths (A2390arcB). 

In Fig. |7| we show the one dimensional rotation curves 
of the galaxies in our sample. These are extracted by sam- 
pling the velocity field with a slit approximately 1.5 kpc wide 
along the major axis cross section, as shown in the recon- 
structed velocity fields of the galaxies (Figs. 2-6). The zero- 
point in the velocity is defined using the center of the galaxy 
in the reconstructed source plane image. The error bars for 
the velocities are derived from the formal 3a uncertainty in 
the velocity arising from Gaussian profile fits to the [On] 
emission in each averaged pixel of the datacube. We note 
that in the one-dimensional rotation curves in Fig. |7| alter- 
nate points show independent data. We also note that the 
data have not been folded about the zero velocity so that 
the degree of symmetry can be assessed in these plots. Hor- 
izontal bars show the effect of seeing, transformed into the 
source-plane. This illustrates how much the velocity field is 
smoothed by the seeing, and how the smearing varies along 
the arc. 

4.2 Comparison to Local Populations 

In Fig |S] we show the TuUy-Fisher relation for the 
high redshift galaxies as compared to pre vious measure- 
ments of the TF in both low redshift JPierce fc TuUvl 
119921: IVerheiienI l200lt iMathewson et al.l Il992l) and high 
redshi ft JVogd Il999t) galax v populations. Whilst IVoeti 
ll999i) and Mathcwson ct al. ( 199 3) us e optic al spectra and 
iPierce fc Tullv. L1991) and .Verheiienl (l200j) use Hi hue 
widths to measure the TF relation, as lMathewson et alJ 
il992h show, with high quality observations the rotation 
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Figure 7. One-dimensional rotation curves from the galaxies in our sample. These are extracted from the two-dimensional source plane 
major axis cross section (and are not corrected for inclination effects). We also show the one-dimensional velocity fields for the two 
galaxies which do not resemble relaxed systems: A2390arcD and A2390arcB. We extract two one-dimensional velocity fields from these 
two galaxies: the first by collapsing the two-dimensional velocity field along the strongest velocity gradient (black) as well as the one- 
dimensional velocity gradient made by collapsing the two-dimensional velocity field along the major axis seen in the [On] emission line 
map (grey). The solid bars represent 0.7" seeing (transformed to the source plane) of each galaxy. 



curve velocity determinations and Hi line widths are tightly 
correlated. Indeed, the scatter between Vopt and vhi is typi- 
cally less than 10 km s"'^ over a ran ge of 60;$t;^300 ( kms~^); 
(Fig. 5 of lMathewson et ahlllflfl^) and therefore it is useful 
to include the local Hi sets in our analysis. 

The four galaxies with well defined rotation curves are 
shown as solid points. As can be seen the points suggest 
a good correlation between magnitude and rotations speed 
similar to that seen in th e local Universe . By a ssuming the 
local slope of the TF from lTullv fc Pierce! i2000l) we measure 
an offset of A/s=0.41±0.34mag of brightening in the rest- 
frame i3-band at fixed circular velocity and place a limit of 
<0.10mag of brightening in the rest-frame 7-band. As we 
discuss below, the small ofi'set in the I-band TF relation is 
consistent with the evolution expected in heirachical galaxy 
formation models. 

In order to plot the remaining two galaxies in Fig|H| we 
must estimate their circular velocity indirectly. A2390arcB 
appears to be a low-mass Hll region, and we place a limit 
of g60kms~^ on any possible velocity gradient (this is also 



comparable to the [On] emission line FWHM). Since we are 
not able to constrain the inclination angle on this galaxy, 
we place this galaxy on the TF relation assuming it is edge 
on, but also include an error-bar which indicates a canonical 
inclination of 60°. In the case of A2390arcD, the galaxy is 
well resolved but does not show a coherent velocity field. In 
this case, we estimate the rotation speed by assuming the 
velocity gradient across the galaxy will eventually develop 
into regular rotation (it is also possible that this velocity 
offsets may develop into a velocity dispersion in an elliptical 
galaxy). The error-bars on the rotation of this galaxy then 
represent the uncertainty as to whether this galaxy is relaxed 
or merging, as well as an inclination correction assuming an 
inclination of 60° (as above). 

This procedure allows us to show that the position of 
these point in the TF relation is reasonable given their opti- 
cal and near-infrared brightness, however, the uncertainties 
in these rotation speed estimates are clearly large. These two 
galaxies should therefore be interpreted with caution. Nev- 
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Figure 8. Left: The arcs from our survey on the TuUy-Fisher relation in rest frame _B-band compared to high redshift sample from Vogt 
et al. (1999). For comparison we show the low redshift local fit from TuUy &; Pierce (2000). The solid points show the galaxies which 
have regular (disk-like) kinematics, but we also indicate the two galaxies for which the dynamics are not well defines as open circles. The 
solid triangle shows field galaxies (L451) at 2; = 1.34 from van Dokkum &; Stanford (2001). Right: The rest frame /-band TuUy-Fisher 
relation compiled from Mathewson, Ford & Buchhorn (1992) (solid points) and from the Ursa-Major Cluster (open circles; Verheijen 
2001). 



ertheless, if these two galaxies are included in the analysis 
we derive an offset of AMs=0.51±0.28 and AM/ <0.10. 

A key gain from studying lensed galaxies is the level 
of detail wfiidi we can extract from their rotation curves. 
For the four galaxies with well defined rotation curves, we 
can also compare the shapes of the rotation c urves with 
those of local (z;$0.03) late-type spiral galaxies bv lCourteaul 
1^93) • This sample is dominated by late-type (Sb-Sc) spi- 
rals and therefore should provide a reasonable comparison 
sample. We fit the observed source-plane 1-D rotation curves 
in our samp le with the ARC TAN and multi-fit parameter 
models from ICourteaul I^Q^)- Using a method the ob- 
served rotation curves are best fit with models with high 
values of 7 and Tt/iopt', (see Fig. |5J. In this formalism, 7 
governs the degree of sharpness of the turn over of the 
rotation curve, rt describes the transition radius between 
the rising and flattening of the curve and Vopt is the ra- 
dius enclosing 83% of the light in the source plane pho- 
tometry. We note that the rotation curve from the arc in 
RGB1745-I-398 is uncertain because the IFU field of view 
only just reaches the turn over, however, the data is suffi- 
cient to place a lower limit on rt /vopt and 7. Such models in- 
dicate slowly rising rotation curves and usually favor fainter 
(or low-surface brightness) galaxies, yet the source frame 
surface brightnesses of the four galaxies are consistent with 
high surface brightness (HSB) galaxies (/it,,o=21.5, 21.2, 20.6 
and 21.0 for A2390arcA, C12236-04arc, A2218arc#289 and 
RGB 1745-|-398arc respectively). In order to investigate how 
sensitive the shape of the rotation curve is to the seeing 
we perform two checks: firstly, we convolve the rotation 
curve with a further 0.7" seeing (transformed to the source 



plane; Fig.|7|l and refit the rotation curve. Whilst the results 
vary between the galaxy rotation curves, the overall effect is 
only modest at best. Typically, we find that the inner rota- 
tion curve flattens the transition radius by Art/ Aropt=0-05- 
0.08, and Ar-t/Ar-opt=0. 05-0.09, and 7=0.18 in the ARCTAN 
and multi-parameter fits respectively. Secondly, we gener- 
ate mock rotation curves for both the ARCTAN and multi- 
parameter fits with low values of rt and 7 and convolve 
these with 0.7" seeing. For the values of rt/ropt and 7 in 
the rotation curves in our data, we find that they could 
be affected by upto Art jr opt =0.06 (arctan), Art jr opt =0.5 
(multi-param) and A7 =0.09. In Figl^jwe therefore show ar- 
row which shows the maximum effect which the seeing could 
have on the measured shapes of the rotation curves. Whilst 
this contributes to the shape parameters we find that in the 
local galaxy population, only about 25% of galaxies have 
similarly shaped rotation curves. This could be an indica- 
tion that the inner parts of the galaxy mass distribution are 
less baryonically dominated than in local spirals and could 
be an indication that the bulge-to-disk mass ratio of these 
high redshift galaxies are lower than those of most present 
day galaxies. 



5 DISCUSSION AND CONCLUSIONS 

In this study we have mapped the two dimensional veloc- 
ity fields of six z ~ 1 gravitationally lensed galaxies using 
IFU spectroscopy. Using detailed mass models for the clus- 
ters we have reconstructed the source morphologies of these 
galaxies. The typical boost in magnitude for these galaxies is 
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z=l arcs 

Local HSB Galaxies 
Local LSB Galaxies 



Figure 9. The distribution of rt/vopt against 7 for the arcs with 
well defined rotation curves in our sample as compared to local 
(z ^0.03) spiral galaxies by Courteau et al. (1997). The arrow in 
the lower right hand corner shows the maximum affect the seeing 
FWHM can have on the shape parameters. The rotation curve 
from the arc in RGB1745+398 is uncertain because the IFU field 
of view only just reaches the turn over and we therefore place 
lower limits on on rt/vopt and 7 for this galaxy. The z=l galax- 
ies in our sample have source-frame brightnesses consistent with 
HSB galaxies, yet the rotation curve shapes are more consistent 
with low-surface brightness galaxies in the local Universe: if the 
maximum effect the seeing can have on the shape parameters is 
taken into account, then only about 25% of local HSB galaxies 
have similarly shaped rotation curves. 



Am = 2, which means the typical (unlensed) i?-band magni- 
tude of our sample of ran ^ 22. Moreover, the gain in spatial 
magnification for a typical galaxy in this paper means that 
0.6" on the sky corresponds to only ~0.5kpc in the source 
frame which allows us to spatially resolve the dynamics of 
these galaxies on much smaller scales than otherwise possi- 
ble. Our IFU spectroscopy of these distant galaxies allows 
us to test whether these high redshift galaxies have regular 
disk-like kinematics, or whether velocity offsets come from 
merging components or other dynamical disturbance. 

With a small sample of high redshift galaxies which 
show bi-symmetric (disk-like) kinematics we can measure 
the evolution of the offset in the TF relation under the as- 
sumption that the slope remains fixed. In Fig. |H| we com- 
pare the rotation velocities and source brightness of the arcs 
which resemble galaxy disks with that of local and other high 
redshift galaxies in rest frame B and /-bands. We also in- 
clude the A2390arcB and A2390arcD in this plot. Whilst 
these two galaxies do not resemble rotating disks and there- 
fore cannot be compared directly to local disks, we show 
the likely range of equivalent circular velocities based on 
the upper limit to the velocity dispersion (A 2390arcB) and 
assuming that the two components are on a merging orbit 
(A2390arcD). These data points are of course highly un- 
certain, not least because we have no way to estimate their 
inclinations on the sky. 

Whilst previous high redshift studies have concentrated 



on the rest-frame _B-band, rest-frame 7-band observations 
provide a more rigorous test of evolution of the TF rela- 
tion since the corrections for dust and on-going star forma- 
tion are much smaller at longer wavelengths iCons elice et alJ 
^005). The rest frame /-band TF therefore gives a clearer 
indication of the true stellar luminosity and hence the ra- 
tio of stellar mass to total halo mass. The position of the 
galaxies on the TF relation in both the B and /-bands 
shows good agreement with local data JPterc^^^Tul^ 
19921: iMathewson et al.lll992l: iHavnes et aljll999l: IVerheiieii 



20011) . This data is in agreem ent with ex isting intermedi- 
ate and high red s hift studies. jlVopt et al.llT 997: Vogt 1999 
Vogt et al. 2002: Zioglcr ot al. 2002; Milvane- Jensen et alJ 



2003c 



2005 



Barden et al.,.,.2 003; ..Bohm et al 2003. : Bamford et alJ 

Con5eUc^^n]2'o05l) . Overall, our observations sug- 



gest a 0.5±0.3mag of brightening in the i3-band TF from 
the local {z — 0) correlation, whilst in the /-band we place 
a limit of < O.lOmag between our z — 1 sample and the local 
z = correlation. 

The theoretical evolution of the /-band TF rela- 
tion from hierarchical models of galaxy formation from 
ICole et al. (.2000i) predict that for any given disk circular 
velocity, the / band luminosity should decrease by ~ 0.1 
magnitudes from z = to z = 1, whilst in the /3-band such 
models predict an increase in luminosity of 0.5 magni- 
tude for the same redshift change. It is useful to compare 
this prediction with the prediction for a simple "classical" 
galaxy formation. We consider a model in which all of the 
galaxy's mass were already in place at 2: = 1, but only half of 
the s tars have yet formed llEggen et al.lll962 : lHopkins et alJ 
l2000h . This toy model produces an evolution in luminosity 
of ~ 0.7 magnitudes from z = to z — 1, during which the 
asymptotic circular velocity remains constant. The small off- 
set of our galaxies from the local /-band TF relation suggests 
a preference for hierarchical rather than the "classical" for- 
mation model. Furthermore, the 0.51±0.28 magnitude evo- 
lution in the rest-frame /3-band is in line with the increased 
star formation activity at z = 1 suggested by stud ies of the 
B-band luminosity function iGiallongo et al.|l2005^ . 

Clearly, it would be dangerous to draw far reaching 
conclusions from a small number of high-redshift galaxies. 
Rather this study should be view as a compa nion to stud- 
ies o f a larger samples of field galaxies (e.g. iFlores et al.l 
'2004*) . The advantage of our lensed study is that we are able 
to resolve rotation curves in detail and clearly identify the 
asymptotic velocity of the rotation curve rather than fitting 
a model rotation curve (based on local galaxies) convolved 
with ground based seeing. In contrast to many ground based 
studies our rotation measurements extend to large radii al- 
lowing asymptotic velocities to be directly measured. 

The present sample is just large enough to place con- 
straints on the evolution of the TF relation, but a signifi- 
cantly larger sample of galaxies are accessible to observation 
using the gravitational telescope technique we have illus- 
trated (there are currently a further ^10 highly magnified 
and spatially resolved z~l giant arcs in the literature which 
are suitable for this study: ISand et a'l]l2005^ . 

However, the principle gain from studying distant galax- 
ies using gravitational telescopes is to study the internal 
properties of the galaxies in detail. In principle, the boost 
from gravitational lensing allows us to achieve an angular 
resolution to galaxies observed with classical techniques at 
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z ~0.1 (e.g. for a typical z ~1 gravitationally lensed galaxy 
in this sample, the spatial sampling of a galaxy is a factor 
of three smaller than unlensed galaxies at the same redshift; 
by contrast this spatial sampling is matched to unlensed 
galaxies at z ~0.15). This increase in resolution allows us 
to study the dynamics, metal abundance and star formation 
of z ~1 galaxies in detail and hence address why the global 
properties of distant galaxies differ from their local counter 
parts. 

In this paper, we restrict our attention to the dynamics 
of the target systems. For those galaxies with well defined 
rotation curves we can use the gravitational magnification 
to compare the shapes of the rotation cu rves with those o f 
local (z;£0.03) spiral galaxies observed bv ICourteaul il997t) . 
We fit the observed source-plane 1-D rotation c urve with 
the AR CTAN and multi-fit parameter models from ICourteaul 
1^93). Using a fit the observed rotation curves are best 
fit with models with high values of 7 and it/iopt- The fits 
from these models show that the high redshift targets have 
slowly rising rising rotation curves compared to typical lo- 
cal galaxies. In the local Universe, such rotation curves are 
associated with fainter, low-surface brightness galaxies, yet 
the source frame surface brightnesses of the four galaxies 
are consistent with high surface brightness (HSB) galaxies. 
In the local galaxy population, only about 25% of galaxies 
have similarly shaped rotation curves and could be an indi- 
cation that the bulge masses of these high redshift galaxies 
is somewhat lower than those at the present day. 

The dynamics of two of the galaxies in our sample 
do not resemble galaxy disks. The first, A2390arcB ap- 
pears to be a highly magnified Hll galaxy, with a veloc- 
ity dispersion of ~ 60kms~^ and a (source frame) ra- 
dius of ~2.5kpc (FWHM). The other galaxy in our sam- 
ple which does not appear to have stable disk kinematics is 
A2390arcD. This galaxy has a disturbed morphology and 
appears to have several components. When combined with 
the velocity field from the [On] emission, we find some evi- 
dence for rotation or interaction. The large [On] line widths 
(;> SOOkms-i FWHM) may originate from either outflowing 
material driven by AGN activity, or as a result of two in- 
teracting galaxies. This second scenario is supported by the 
fact that we also identify a bright knot, or companion offset 
by 5kpc and 480±60kms~^ in projection; the most likely 
interpretation of the system is that this galaxy has recently 
undergone a tidal interaction/merger which has produced 
the disturbed and complex morphology. The selection cri- 
teria for our analysis was simply that the galaxies must be 
highly amplified, and therefore it is interesting to note that 
although we only have a small sample, of the six galaxies 
we have studied, two do not have disk-like kinematics or 
morphologies, giving some indication as to the mix of stable 
disk- like galaxies to currently assembling galaxies at z 

The next step in this study is to complement these ob- 
servations with similar observations in the near-infrared. At 
z H/3, [Olll] and Ha are redshifted to beyond l^m and 
until recently have been inaccessible to Integral Field Spec- 
troscopy (IFS). However, recent developments in designing 
efficient image slicing IFU's allow us to probe beyond l^m, 
and therefore we can combine the rest-frame optical emission 
lines to probe the distribution of reddening t hrough the well 
studied R23 index; 7l23=([Oll]-|-[Olll])/H/3 JZaritskv et alJ 
I1994I) , as well as the reddening corrected star formation rates 



through the Ha:H/3 decrement. For those galaxies which we 
have identified as disk-like at z ~1, such observations will 
give unique insight into the processes of disk galaxy forma- 
tion, and give important constraints on the assembly on the 
disks of present day galaxies. 

The gain in spatial resolution which we have achieved 
by using a gravitational lens to magnify and stretch distant 
galaxies demonstrates the science that will soon be possible 
with Adaptive Optics Integral Field Spectroscopy (AO-IFS) 
on (non-lensed) galaxies ai z — 1 (e.g. OSIRIS on Keck, 
NIFS on Gemini and SINFONI on the VLT). These stud- 
ies will probe the structure of young galaxies on lOOmilli- 
arcsecond (mas) scales, although long exposures will be 
required to compensate for the lack of lensing amplifica- 
tion. For a non-lensed source ai z — 1, this corresponds 
to 600pc. However, the greatest gains will come from com- 
bining AO-IFS with the powerful gravitational lensing tech- 
nique. For a lensed galaxy with magnification factor of ten 
in flux (a factor of three in linear scale), the effective diffrac- 
tion limit of the telescope is reduced by a factor three and 
lOOmas corresponds to only 200pc in the source frame of the 
galaxy, sufficient to resolve the star-formation, kinematic 
and chemical properties of individual Hll regions (which 
typ ically have sizes less than 400pc in local spiral galax- 
ies; iGonzalez Deleado fc Pere'^ll9^7^ . We look forward to 
undertaking these studies in the coming years. 
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